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THE  ASYMPTOTIC  VARIANCE  OF  HIGHER  ORDER  CROSSINGS 


By  Benjamin  Kedem  and  George  Reed 


1.  Introduction 

The  computation  of  the  variance  of  the  number  of  axis  crossings 
by  a  stationary  process  entails  technical  difficulties  in  both  the 
continuous  and  the  discrete  time  cases.  What  makes  this  problem 
so  intractable  is  that  the  variance  is  expressed  via  quantities  which 
involve  conditional  intensities  and  fourth  order  orthant  probabilities 
which  are  difficult  to  compute.  This  is  so  even  in  the  Gaussian 
case.  See  Cramer  and  Leadbetter  (1967),  Sec.  10.7,  Cox  and  Isham 
(1980)  pp. 33,87,  Kedem  (1980)  p.68.  On  the  other  hand  when  the 
process  is  sufficiently  oscillatory,  it  is  possible  to  estimate  the 
variance  under  suitable  conditions  by  a  simple  expression.  This  is 
shown  here  by  examining  the  variance  of  the  higher  order  crossings. 

The  notion  of  higher  order  crossings  has  been  introduced  by 
Kedem  and  Slud  (1981),  (1982)  as  follows.  Let  }~__00  be  a 

stationary  discrete-time  process  and  define  a  sequence  of  binary 


processes  by 


1  if  Vk"1Zt  >  0 
0  otherwise 


k  =  1,2,..., 


where  V  is  the  backward-difference  operator,  VZ1 
Then 

„  _  N;1  ,r„(k)  ,,  „<k), 

(2)  DkjN  =  lCXt+1  *  Xt  ] 


Zt  -  Zt-1 


2 


is  the  number  of  axis-crossings  by  { >  and  is  called 
the  (number  of)  higher  order  crossings  of  order  k  associated  with 
{Z^}.  The  main  result  of  higher  order  crossings  is  (Kedem  and 
Slud  (1982))  The 

Higher  Order  Crossings  Theorem  (HOCT):  If  {Z^}  is  a  strictly 
stationary  process  with  finite  variance  and  P  (Z^  5  Z^tft)  =  0  and 
if  the  spectral  distribution  function  F  satisfies  the  condition 
dF(ir)  >  0  then 

f,  v  (...01010101...  with  probability  1/2 

a)  <x;k;)  =»< 

T  1  ...  10101010 .. .  with  probability  1/2 

(ii)  lim  lim  N-1D,  M  =  1  with  probability  1. 

k+oo  N->«»  K » N 

Thus  we  see  that  N  eventually  increases  with  k  for  sufficiently 
large  N.  This  monotone  property  of  higher  order  crossings  can  be 
associated  with  the  amount  of  information  they  carry  as  exemplified 
in  Kedem  and  Slud  (1981),  where  a  discrimination  statistic  made  of 
higher  order  crossings  is  introduced.  However  the  HOCT  does  not 
tell  us  the  rate  of  increase  in  the  ^ .  The  present  paper 

provides  a  refinement  of  the  HOCT  in  that  we  establish  under  some 
conditions  the  rate  of  decrease  in  Var(D^  ^)  as  k  -*■  00  for 
fixed  N.  The  fact  that  this  variance  should  decrease  is  apparent 
from  the  HOCT  since  the  probability  of  two  consecutive  l's  in 
(X^)^  approaches  zero.  More  precisely,  define  the  intensity 
associated  with  {X^.k^} 


3 


Then  our  main  result  states  that 


Theorem  1.  Let  {Z^}  be  an  m-dependent  stationary  Gaussian 
sequence.  Then  for  fixed  N 


(4) 


lim 

]<-*■« 


var(iy; 

<N-m^k) 


i. 


The  limit  (4)  certainly  holds  for  any  finite  moving  average.  An 

explanation  of  (4)  is  provided  by  the  intuition  that  by  the  HOCT, 

{X^.  '  }  "becomes  Markovian"  and  A^  -+•  0  as  k  -*  °°.  Note  that  if 

{X^k^}  is  truly  Markovian  then  VarCD^  ^)  =  (N-l) A^k^ ( 1-A^k^ ) 

(k)  * 

which  yields  (4)  when  A-^  approaches  zero.  See  Kedem  (1980) 
Chapter  3 . 

k-1 

Since  for  moderate  and  large  k,  {V  Z becomes  oscillatory, 
i.e.,  Dk  n  lar2e>  apparent  from  the  HOCT  that  (4)  sheds 

a  great  deal  of  light  on  the  behavior  of  the  variance  of  N 

(number  of  axis  crossings)  for  oscillatory  processes.  It  should 
be  noted  that  (4)  is  not  entirely  intuitive  since  in  general  we  have 


Var(Dk,N) 

<N.mR 


(N-l)Var(- 


D 


k,n, 


N-l' 


TkT 


(k) 

(N-l)(l-A^  ) 


for  all  k,  and  for  N  >  2.  Other  than  N  =  2,  this  general 
bound  has  little  resemblance  to  (4).  On  the  other  extreme,  many 
simulations  with  low  order  autoregressive  moving  average  processes, 
at  least  for  small  k,  suggested  that  the  limit  should  be  close  to 
zero.  It  is  only  as  k  increases  that  (4)  holds.  It  seems  to  us 
therefore  that  for  large  k,  in  light  ot  our  earlier  comment,  the 


assumption  that  { X^.  }  is  approximately  Markovian  is  not  entirely 

unrealistic  as  it  immediately  leads  to  (4).  In  this  paper  however 
this  assumption  is  not  made. 

The  proof  of  Theorem  1  is  carried  out  by  first  proving  a  series 
of  lemmas.  In  Section  2  we  establish  a  two  sided  inequality  while 
in  Section  3  we  construct  certain  continuous  differentiable  functions 
which  establish  the  rates  of  convergence  of  the  correlations  of 


5 


2 .  An  Inequality 

In  order  to  prove  and  motivate  Theorem  1  we  need  to  establish 

a  few  basic  facts  concerning  axis  crossings  by  Gaussian  sequences. 

( k ) 

From  now  on  it  will  be  convenient  to  suppress  k  in  A^  and  in 
(k) 

and  we  shall  simply  write  A  ^  and  Xt  with  the  understanding 

that  both  quantities  depend  on  k. 

First  observe  that  D.  ..  can  be  written  as  a  sum 

k  ,N 


d2  *  d3  +...♦  dN 


where 


(6)  dt  =  Xt  ♦  Xt.1  -  2XtXt_1 

and  is  also  a  function  of  k  since  X^  is. 

Lemma  1.  Let  {Z. }”  be  a  stationary  Gaussian  process.  Then 

—  t  L - 00 

uniformly  in  t 


Cov(Xt  j>dt)  =  °>  D  =  0,±1,.. 


Proof.  The  Gaussian  assumption  implies  that 


Cov(Xt..,Xt+Xt.1-2XtXt.1)  =  [|  X.-i]  *  [i 

“  2[i  Ax  +  i  A.  +  £  \._±  -  =  0 

by  noting  that  A.  =  |  ^  Sin"1  p  .  where  p.  is  the  correlation 

function  of  {Zt  > .  (See  Kedem  (1980),  p.33.)  The  orthogonality 
eolation  (7)  implies  that  and  are  independent  since 

both  variables  are  binary.  We  can  therefore  rewrite  Lemma  1  as 


Lemma  1  * .  Let  be  a  stationary  Gaussian  process.  Then 

(8)  Cov(g(Xt_;. )  ,dt)  =  0,  j  =  0,±1,... 


for  any  function  g. 

Next  we  examine  the  relation  between  X.  .  X^  and  d^. 

t-D  t-3+1  t 

Unlike  the  orthogonality  relation  (7)  the  covariance  between  the 
last  two  quantities  cannot  be  evaluated  directly  but  may  be 
approximated.  Without  loss  of  generality  we  will  focus  on 
CovCd2,XtXt_1). 

Lemma  2 .  Let  {Z^}”,  ^  be  a  stationary  Gaussian  process.  Then 
we  have  the  two  sided  inequality 


/n\  ,  1  .  1  .  1  ,  .  1  ,2  1.  I,  1.1,2. 

(9)  max(^  i  _  4  ^t-3  "  2  Ai  +  2*  ^l’~4  At-1  +  4"  At-3  “  2  A1  +  2  Ai' 


S  Cov(d2,XtXt^)  5  ±  A*. 


Proof.  First  note 


Cov(d2  >XtXt_1) 


Fd2XtXt-l  '  Ed2EVt-l  5  (EXtXt-lH1-Ed2) 


1 

2 


X 


2 

1  * 


Fcr  the  left  hand  side  we  have  four  cases  corresponding  to  the 
"removal"  of  either  X±  or  X2  or  Xt  or  X^.  Thus  for 
example  if  we  "remove"  X^  we  have 

-2EX1X2Xt  xXt  >  -2EX2Xt_iXt  =  At-3+^  At-2  +  J  X 


from  which  follows  that 


Cov^-Vt-l’  ?  +  f  Xt-2J  *  At-3] 

+  [|-|  Xt-3'f  Xt-2'f  Xl*f  Xl] 

*  I  Vl-iVj  -  f1!*?1! 


and  similarly  for  the  other  three  cases. 


Now  observe  that  in  general 


(10) 


Var( ) 


x1(i-x1) 


and  therefore  by  Lemma  1,  Theorem  1  is  proved  if  we  can  show  that 
Cov(d  ,X  X  )  =  Cov(d  )d  )  =  o(X..)  for  arbitrary  t,s.  However 

S  X  X  S  t  X 

this  is  not  immediate  from  Lemma  2  and  requires  some  careful  con¬ 
siderations  to  be  discussed  rigorously  in  the  following  Section. 

It  should  be  made  clear  though  that  (9)  is  instrumental  in  proving 
that  Cov(d  ,d  )  =  o(X  ).  We  can  support  this  claim  by  resorting 

S  L  X 

to  an  example  heuristic  in  nature. 


Example  1 . 


and  so 


T  1  * 

Assume  X^  is  Markovian.  Then  X^  =  y  +  ^ 2X^-1) 3 


4VAt-l  t-3^ 


\  X1(X1-l)(2X1-l)t‘3, 


Then  the  left  hand  side  of  (9)  becomes 


Cov(d2,XtXt  l)  >  max(|  X-^ X^-l) ( 2X^1) 


t-3.  I  x  +  1  x2 

2  A1  2  Al* 


-ix,(x,-l)(2x .-l)1"3  -  i  X  +  i  \h  =  o(X,  ) 


since  — 2  X^  is  always  eliminated  from  one  of  these  las'*-  two 
expressions  regardless  of  whether  t  is  odd  or  even.  Thus  for 
a  Markovian  X  (9)  and  Lemma  1  imply 


(11) 


Cov(d2,dt> 


o<V 


Another  indication  that  (11)  is  plausible  is  furnished  by  the  next 
numerical  example. 


Example  2 .  Assume  {Zt}  is  a  Gaussian  white  noise.  Then 


(12) 

and 

(13) 


f  *  ?  sin'1,-m) 


W(k7I)2,J( 


2k+l 

(k+1) 2(k+2 


;>)■ 


Incidently,  here  we  see  the  direct  dependence  of  X^,X2  on  k. 
Now  define  X*  by 


(14) 


We  maintain  that  as  k  -*■  ®  X*  approximates  X2  very  closely. 
This  can  be  best  seen  from  the  following  calculations. 


X2 

X* 

A2 

h 

0.97155556 

0.97154135 

0.001422932 

0.98726994 

0.98726865 

0.00636567 

0.99099766 

0.99099720 

0.00450140 

0.99715295 

0.99715295 

0.00142353 

0.99909968 

0.99909968 

0.00045017 

X  -A* 


0.00001421 


0.00000129 


0.00000046 


0 


0 


We  see  that  as  k  ■*  00  the  white  noise  case  implies 


(15) 


A*  +  o(A1) 


This  can  also  be  verified  directly  from  (12)— (13) .  It  follows 
that  as  k  -+  00 


|  -  Sin_1(-l+2p.2)  +  o(\,) 

L  IT  11 


(16) 


1  +  —  Sin  1(-cos  2tt(A  -4))  +  o(A  ) 

2  it  12  1 


7  -  +  o(V 


1  -  2A1  +  o( A1) . 


But  in  general,  for  any  Gaussian  stationary  sequence 
(17) 


Cov(dt>dt+1) 


and  so,  in  the  white  noise  case 

Cov(dt,d^+1)  =  oCA^)  -A2  =  ©(A^). 

This  example  indicates  that  we  may  also  have  Cov(dt,ds)  =  o(A 
s  i  t,  for  any  two  terms  in  (5).  However  the  proof  of  this 


assertion  requires  the  establishment  of  some  further  facts  and 
the  use  of  (9). 


3.  Proof  that  Cov(d. ,d  )  =  o(>.) 

_ t  s _ i_ 

As  stated  in  Section  2  it  is  sufficient  to  show 
Cov(d2>XtXt  =  o(A^)  in  order-  to  complete  the  proof  of  the 
theorem.  From  (9)  it  is  clear  that  CovCdjjX^X^  )  «  oCA^)  go 
the  problem  reduces  to  showing  o(A^)  s  Cov(d2 ,X^ , X  )  using 
the  left  hand  side  of  (9). 

(k ) 

It  is  more  convenient  at  this  point  to  express  A  in  terms 

of  the  correlations  of  the  k^*1  differenced  process,  From 

Kedem  (1980)  we  have  A^k^  =  i  +  —  sin  1p.  Ignoring  the  A^ 

t  2ir  t  ,k  &  b  1 

term  in  (9)  and  suppressing  the  k,  it  will  be  sufficient  to  show 


(18) 


.  .  ,  1  -  -I  1  •  -1 

lim(w  sin  p 47  sin  P 

k-*-°° 


t-3 


1  1  .  -1  ,,,1.1  .  -1  V 

4‘Hsln  pl,/(2+isin  01> 


=  0 


or 

(19) 

,  .  ,  1  .  -1  .  1 
lim(-—  sun  P-t-].  +  4tt” 
k-*-°° 


.  -1 
sm  p 


t-3 


1 

4 


2n 


.  -1  W,1  .  1  .  -1  V 

sxn  p1)/(— +-  sin  p.  ) 
1  2  it  1 


=  0 


The  solution  involves  two  ideas.  First  treating  all  functions 

of  interest  as  functions  of  p  ,  since  this  is  the  basic  variable 

1 ,  K 

at  hand  as  is  apparent  from  the  HOCT,  and  secondly,  circumventing 
the  problem  that  both  numerator  and  denominator  in  (17),  (18)  go 
to  zero  by  constructing  an  equivalent  limit  with  continuous  differ¬ 
entiable  functions  and  using  l'Hopital's  rule. 

The  following  lemmas  will  show  how  to  construct  the  equivalent 

limit  by  replacing  p  .  with  an  appropriate  function.  First 

t  *  K 

define  the  following  notation  associated  with  the  k'th  and 
(k+D'st  differences. 


7,k+l  KV,k 


V,k 


Pl,k+1-Pl,k 


=  C^‘[-l,p1  ,]  =  set  of  continuously  differentiable  functions 
1  *  1 


on  the  interval  [-l,p.  ,3 

1  j  1 

Again,  as  above,  whenever  it  is  not  confusing  the  subscript  k 
is  suppressed. 


Lemma  3.  If  lim  Rn  ,  exists  then  for  each  sequence  of  ordered 
k-°°  V,k 

pairs  {(p,  ,,p  .)}*  ,  r  fixed  integer  >2,  there  exists  a 

Ijk  TjK  K“1 

function  p  ( p, )  €  C1  such  that 
r  1 


(i)  pr(pi,kJ 


(ii) 


J  i  • 

r,k 


dp. 


limA  =  a  is  independent  of  the  function 

01~lt  dPl 


constructed. 

(iii)  a  satisfies  the  difference  equation: 


a  ..  -  2a  ( l-2a.. )  +  a 
r+1  r  1  r-1 

,  •  dpV 

a*  =  Pl-l+  dpl' 


( — 1 ) 1 2  where 


n. 


Note  that  p,  is  p,  .  but  at  the  same  time  it  is  used  as  the 

X  1  j  K 

variable  of  differentiation.  Also,  a1  is  defined  differently 
from  a0 ,a  , . . .  • 


from  the  pairs 


Proof.  We  first  construct 


co 

{( Pf  k»Py  and  then  the  °"ther  functions  will  be  produced 

inductively. 


By  a  direct  computation  we  have  lim  p  =  -1  and  it  has  also 

fc->-oo  1»* 

00 

been  shown  that  {p.  ,  ,  is  a  strictly  decreasing  sequence 

x  j  K  K“ 1 


(Kedem  (1982))  to  -1.  In  what  follows  the  two  equivalent  notions 


of  k  -*•  00  and  -*•  -1  are  used  interchangebly  without  loss  of 


generality. 

QO 

Since  (p,  is  strictly  decreasing,  for  any  finite  set  of 

J.  5  K  K  —  X 


values  { ( p  ,p„  )}[J  we  can  construct  a  function 

*  j  K  K- X 


n*n)(p_)  €  C1[p1  ,p.  ,]  such  that 

v  l  i,n  i,i 


nvn,<‘>i,k) 


'V,k 


and 


dh 


(n) 


dp 


(p  )  =  R 


l,k 


If  we  let  nyI1+1^(P1)  be  the  extension  of  n^^Cp^  then  we  can 


define 


and 


w  = 

lim 

n->-°° 

( n ) .  v 

^  1  ^ 

for  Px  6 

p7(-l)  = 

lim 

k-*“ 

PV(pl,k5 

=  lim  +  py(p. ) 

P-j^-1 

dP  n 

^-X(-l)  = 
dpx 

lim 

k-*-<» 

RV,k 

lim  ,  -r— 
d°i 

By  the  construction  of  n^n\p  ),  p_(p  )  e  C'(-l,p  ]  and 

v  l  vi  .  l,i 


since  p„(-l)  is  well  defined  by  the  HOCT  and 


dP< 


dp 


-(-1)  is  well 


defined  by  our  assumption  lim  ,  exists,  then  p„(p,)  €  C 
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dP, 


To  see  that  lim  .  is  independent  of  the  construction 

P-^-l  dPl 


of  p  (p  ),  note  that  any  function  satisfying  (i)  must 

do7 

satisfy  lim+  —  =  lim  R^  .  .  This  is  so  simply  since  by  con- 

p-*— 1  P1  k-*-°°  *  , 

dpr 7 

struction  for  every  e  3K  ?  for  k  >  K  I  R„  ,  - -j— -(p„  ,  )|  <  e. 

1  v,k  dp^  l,k  1 

And  so  lim  R0  ,  is  independent  of  the  construction  of  pn. 

k-K»  v » K  v 

We  can  now  define  the  functions  p  (p,).  Since 

r  1 

Py  =  (2p1-p2~l)/2( 1-p^)  we  have 


*2  ,k  =  "2(1_pl,k)pV,k  +  2pl,k  ‘  li 


So  we  can  define 


(20) 


P2(pi)  =  -2(l-P1>Pv(p1>  +  2px  -  1 


In  general,  assume  the  existence  of  pr  ^(p.^)  jp^^p^)  SC1.  Then 


from  the  fact  that  p_  =  (2p  -p  ,-p  , )/2(l-p,)  we  obtain 

V  r  r+1  r-1  1 


r,k 


'2<1_pl,k)pV,k  +  2pr-l,k  "  pr-2,k 


and  we  can  define 


(21)  p^p^ 


.r-1, 


-2(l-p^)p^  (p^)  +  2pr_^P-j^  ”  Pr-2^P1^ 


where  p^r~^(p^)  €  can  be  defined  since  we  can  always  construct 


a  function  p^  tp^)  using  the  same  sequence  of  pairs  and  the  same 

( r*  1 ) 

constructions  that  we  used  for  p^  ^(p^),  so  ^et  py  (p^)  = 

(p^( p^)  )  a  composition  of  C  ^  functions. 


Hence  for  each  r  there  exists  p^(  o,  )  satisfying  (i)  and 
(ii).  To  generate  the  difference  equation  in  (iii)  consider 


pr*l( °1> 


Differentiate 
dpr+l  . 


dp, 


(r ) ,  x  .  \  “"V  .  0 


dp 


(r) 


2p— (Pl)  -  2(1-Pl)  dpi 


dP. 


+  2 


r-1 


dp,  dp. 


and  take  the  limit  as  p^  •*  -1 


C  22) 


r+1 


2 C  —  1 )  -  4 (  lim 


dp 


(r) 


.  ,  -)  +  2a  -  a  ,  . 

,  +  dp,  r  r-1 

p,-*--!  1 


dp 


(r) 


We  can  write 
(r) 


dp. 


<F>  dPy 
dp,  dPa 


but  note  that 


fp; 

+  dr 


dP, 


dp 


(r) 


lim  .  =  lim  ,  or  lim  .  — ~ —  =  a  a,  and  therefore 


Pl-i  pv 

(22)  becomes 


Pl-i-i 


+  dp. 


Pt  "*'_1 


V_ 

+  dp. 


r  1 


r+1 


2 ( — 1 )  -  4a  a  +  2a  -  a 

r  1  r  r-1 


which  can  be  rewritten  to  give  the  difference  equation  in  (iii). 
This  completes  the  proof  of  Lemma  3 . 

r+i  2 

Lemma  4.  If  lim  Ry  ,  =  1  then  af  =  (-1)  r  . 

k-*00  * 


Proof.  lim  Ry  k  =  1  means  a 

k->®  » 

becomes : 


1  so  that  the  difference  equation 


ar+l  + 


2a  +  a  . 
r  r-1 


(-l)r2. 


A  specific  solution  is  a 


r+1  2 

(-1)  r  .  The  general  solution 


to  the  homogeneous  equation  a  ,  +  2a  +  a  =  o  i< 

r-±  r  r-1 


ot(-l)  +  8(-l)r  r* 

But  using  (20)  a2  =  -4  and  using  (21)(or  simply  (22)  with 

r  =  2)  with  r  =  3  a^  =  9.  These  conditions  result  in  a  =  8  =  0 

and  thus  a  =  (-l)r+^  rM 
r 

Lemma  5 ♦  lim  R„  =  1  if  the  process  is  m-dependent. 
k-*°°  »K 

Proof .  The  result  can  be  derived  from  the  expression  for  p 

1 5  k 

given  in  Kedem  and  Slud  (1981): 


(23)  p. 


-CKD-ai-- 

Pm — PH - 

(  I-  2p  [  )+...+  (-1) 

\  k  \k-l 


+  (-1)  P, 


For  white  noise,  it  is  a  straightforward  computation  since 

P1  k  =  ”  k+T*  In  2eneral  Pi  =  0  for  i  >  m,  let  k  >  m  so 

that  there  are  no  new  correlation  terms  in  (23)  as  k  increases. 

Then  p^  -  p^  ^  can  always  be  written  as  the  quotient  of  two 

finite  polynomials  in  k:  P1(k)/P  (k).  Since  there  are  no  new 
terms  in  the  expression,  as  k  increases,  then  p  .  -  p 

V  )  JC*  X  V  j  K 

will  be  the  quotient  of  the  same  two  polynomials  only  in  k  +  1: 

P.(k  +  1)  P  (k) 

P1(k+1)/P2(k+1) .  Hence  ^  =  p-  —  •  p~^+yy  and  since  these 

are  finite  polynomials  in  k,  lim  R„  ,  =  1 

k-*°°  ’k 

With  these  results,  it  is  now  possible  to  evaluate  (18)  and 
(19).  If  t  is  even  consider  (18).  Using  Lemma  3  we  can  write 
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'*  m 

c: 


<i 


i 


c 


i 


i 


i.y 


an  equivalent  limit  by  replacing  pt  R  with  pt(p  )  and  considering 
(24) 


lim1+(4TT  sin  pt-l<pl)  '  411I  Sin 
Pi  **■“•*• 


1  1  .  ,  -1  s  ,,1  A  1  .  -1 


pt-3<pl>  •  4“  27  Sin_  pl)/(2  +  ^sin* 


Using  l'Hopital’s  rule: 

1  1 


(25)  lim. 


A-“li 


dpt-i  l 
"  4 


dp 


dp. 


L  2  dpl 
/1_Pt-3 


t-3  1  _ 1 

2 


px*-i 


/C 


This  limit  can  be  evaluated  if  we  assume  m-dependence.  Using 

2 


P1  -  ^p  1  1 
Lemmas  4  and  5  we  note  that  lim  +  - ^  =  lim  + 


Hence,  lim  + 
>1 


l  2 

/  1~pl  _  1 


n  -*•- 1'  1-a^  o  -v-1 '  dpt  Ta?T' 

P1  Pt  P1  -2p  1 

t  dpl 


.  Thus  (25)  becomes: 


/Iati 


(26) 


i  3t-l  1  at-3  1 

“/rvir  4/^'7' 


When  t  is  even 


it_i  =  +<t-l)%  at_. 


+  ( t- 3  ) 


and  (26)  equals 


i(t-l)  -  jj(t-3>  -  i  =  0. 


When  t  is  odd  a  similar  computation  shows  (19)  reduces  to  zero 
as  well.  So  for  all  t  Cov(d2 =  o(\^).  This  completes 


the  proof  of  the  theorem. 


£.my. ^  _■« _* , —  __  *, " ,  . , 
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As  the  m-dependence  assumption  was  used  only  in  showing  th.it 
,  -*■  1,  k  -*•  it  is  clear  that  the  theorem  can  be  extended 

to  cover  a  more  general  case  for  which  this  requirement  still  holds 
provided  approaches  zero  fast  enough.  We  conjecture  that 

Theorem  1  holds  true  at  least  for  stationary  autoregressive  moving 
average  Gaussian  sequences.  The  proof  of  this  requires  a  more 
careful  examination  of  (23). 
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